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Device, method and system for measuring the distribution of 
selected properties in a material 

Field of invention 

This invention relates to a device for measuring the 
5 distribution of selected properties in a material, and in 

particular a device that non- contacting and non-destructively 
measures the spatial distribution of material properties, such 
as density, water contents and temperature of materials, by 
detecting electromagnetic radiation. The invention also 
10 relates to a method and a system. 

Background of the Invention 

Many industrial processes depend on the measurement of 
material properties as temperature, water contents and 
material density. A close monitoring of these material 
15 properties results often in increased efficiency and improved 
product quality. Additional benefits are likely to occur, if 
such measurements can be accomplished fast and in a non- 
destructive, non -invasive, and non- contacting way with 
acceptable accuracy. 

20 As an example the determination of the temperature 

distribution in foodstuff during heating process. Here, an on- 
line monitoring of the temperature distribution helps to avoid 
cold spots where bacteria are not eliminated completely or to 
reduce the overdue heating time spent to ensure complete 

25 bacteria elimination. This results in a reduced heating time 
and reduced energy consumption as well as in an increased 
throughput of the production line. 
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Material properties are traditionally measured by some form of 
destruction (sample separation, peeking ) but can often be 
measured by the analysis of transmitted electromagnetic 
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radiation by evaluating the dielectric response of the 
material. Measurements using electromagnetic radiation are 
generally contact-free and non-destructive. 

A suitable frequency region of electromagnetic radiation to 
5 determine material properties as temperature distribution , 

water contents and density is the lower microwave region where 
water absorption is not too large and the wavelength is 
already short enough to ensure reasonable spatial resolution. 
The determination of the above material properties is achieved 

10 by analysing the dielectric response of the material based on 
the material's polarisability . Dielectric data of a material 
sample are typically obtained in analysing the electromagnetic 
wave's reflection and transmission properties or a combination 
of both. In order to obtain a distribution of the material 

15 properties, a three-dimensional image of the material's 

dielectric response must be measured. This requires to move 
the microwave detector setup and the material sample relative 
to each other. 

Prior art instruments make use of either a single measurement 
20 frequency or the emission frequency is swept within a 

frequency interval (FMCW) and the average delay time is 
calculated from the obtained data. 

Prior art that approaches to dielectric imaging, use the 
transmission of electromagnetic radiation of a single 

25 frequency (or a small band) between a multitude of antenna 
locations or the material sample is shifted and rotated or 
shifted in two dimensions in order to obtain a spatial 
resolution. Based on these data the dielectric image is 
obtained e.g. by the well-known CSI (contrast source 

3 0 iteration) method where the location and strength of 

polarisation sources are obtained in an iterative process . 
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Using techniques well-known to a skilled person (i.e. Contrast 
Source Iteration CSI, as described below) the electromagnetic 
picture is used to calculate the unknown dielectric functions 
in the dielectric picture. 

5 Starting from Maxwell's Equations (as described by R.F. 
Harrington, in the book with the title w Time Harmonic 
Electromagnetic Fields", published by Mc . Graw Hill 1961) one 
assumes that any region where the dielectric function is 
different from unity, the electromagnetic field creates bound 
10 charges due to polarisation. These bound charges are created 
by the electric field itself and they oscillate with it 
resulting in an additional current component: 

j(p)= s{p)u{p) 

15 where the current density is j, the electric field is u and 
the dielectric function of the material is s and of the 
background is denoted e b . Assume p and q to be two position 
vectors in a two dimensional cross section of the measurement 
gap. D is a domain which contains the cross section of the 

2 0 material sample. The vector q denotes the source point of the 
electromagnetic radiation. Based on that a general relation 
for the connection between the electric fields in the 
measurement space is obtained formally by applying the 
definition of a Green's function for the electric current: 

j 

2 5 Uj(p) = k 2 \G{p,q) j{q) dv{q) 

D 

Inserting the above current density relation and splitting the 
integral yields : 
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",(/>)= * 2 \G{p y q) s b u(p) dv(q) + k 2 \G{p iq ) [e{py€^ u{p) dv(q) 

D D 
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Here the first term denotes the electric field when the 
dielectric response of the background is present only, the 
second term stands for the fields generated by polarisation 
i.e. a dielectric contrast. The fields when only a background 

5 is presented are referred to as incident fields u™ . Then the 
field at an observation point incident from the radiation 
source is (according to an article by P.M. van den Berg, B.J. 
Koo j , R.E. Kleinman, with the title "Image Reconstruction from 
Iswich-Data III", published in IEEE Antenna and Propagation 
10 Magazine, Vol.41 No. 2 April 1999, p. 27-32) : 

Uj(p) = < 00 + * 2 l G (P><!> X{q) u s {q) dv{q) 

D 

where G denotes the two-dimensional Green's function of the 
electromagnetic problem 

15 G{p, q )=^(k\p- q \) 

and the polarisability function x depends on the dielectric 
function of the material e and the background e b in the 
following way: 

Defining scattered fields f one obtains directly: 
F,(r)=« > (r)-«f(r)=A I \G{r,q). X {q) ■«,(*)• Ml) 

D 

25 From this an integral equation for the scattered electric 
field at any point r is set up. 

F,.(r) = -U J ^)h( 9 ) + «f( g )]rfv( ? ) 

^ D 
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This relation is fulfilled exactly when r is equal to the 
antenna location and the F t (r) are measured values of the 
scattered fields for a given wave vector k for a frequency f : 
In 

k = /. The values of F,(r) for the points interior to the 

c 

5 region D are only fulfilled approximately. So the above 

relation has to be solved for a set K of k vectors and a set Q 
of internal points resulting in a [^*£?]x[A^0] non-linear matrix 
problem for the fields F,(r) and the polarisabilities x( r ) • 

In matrix form the state equation becomes: 

10 u^U™ + G%U 

whereas the frequency relation is: 
F = Gxu 

15 Introducing the contrast sourcep = ^-« the above relations 
become 

<p = xu inc + xG<p at all Q interior points, for any of the K 
measurement frequencies 

F = G<p at a single antenna location, for any of the K 

20 measurement frequencies 

Using the method of conjugated gradients sequences for the 
contrasts and the contrast sources solving the above problem 
are obtained. 

Objective of the Invention 

25 A device has been designed to measure the spatial distribution 
of the temperature, water contents and density distribution in 
a material based on the dielectric and magnetic information 



contained in transmission measurements obtained using 
microwave radiation. 

This invention covers two methods to resolve such information 
from measured data: 

(a) The temperature, density and water contents profile can be 
obtained by interpolation between a set of previously measured 
material samples where the profiles are known in advance. 
There the measurement result is found by a best fit to the 
interpolation database. 

(b) The said profile is found by direct calculation of the 
inverse scattering problem resulting in a known distribution 
of the dielectric and magnetic properties. Based on models on 
the dependence of the dielectric and magnetic properties as 
functions of the wanted parameters, a map of the said 
properties is obtained directly. 

The instrument proposed here may only use one mechanical 
scanning dimension. Due to the usage of a multi-channel 
antenna and a multitude of frequencies, a two-dimensional 
cross -section of the dielectric picture is obtained . This 
calculation process involves a novel method related to 
contrast source iteration where the location and strength of 
polarisation sources are obtained in an iterative process 
based on transmitted electromagnetic field measurements at a 
multitude of frequencies. Thereby the antenna patterns must be 
frequency dependent and they are assumed to be directed in 
cross section of the sample allowing an essentially two 
dimensional approach . 

In order to facilitate the calculation of the dielectric 
parameters, regions where the dielectric properties are at 
first order constant are obtained by an e.g. evaluating video 



pictures taken from at least two different points of view with 
overlapping image region. From these video pictures a 
reasonable guess of the material sample's dielectric structure 
is made. As an alternative ultrasound images can be used for 
the same purpose or a three dimensional image of the material 
sample may be stored in a memory. 

Summary of the Invention 

The object with the invention is thus to provide a device that 
measure the spatial property distribution in a non- contacting 
and non-destructively way. 

In accordance with the invention this object is achieved by 
the features in the characterising portion of claim 1 . 

This object is also achieved by the features in the 
characterising portion of claim 12 and claim 17, 

An advantage of the invention is that it provides on-line fast 
measurements of spatially resolved material parameter 
distributions by means of a combined application of microwave 
reflection and transmission measurements and a three 
dimensional contour of the material . 

Traditionally the temperature and density of the material 
samples is obtained by probing a certain fraction of the 
material samples. This method allows a complete on-line 
monitoring of all material samples in production increasing 
the degree of product control. 

The accuracy of the measurement is checked by means of 
calibration samples with known constituents and known 
temperature profile which are measured at regular intervals. 
Thereby it is sufficient to perform invasive temperature 
measurements after the sample has been measured at different 
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points of the sample and compare them to the instrument's 
findings. As a additional verification process the same 
procedure can be repeated when the sample has e.g. cooled 
down. 

5 Summarising the advantages of this invention it provides a 
non-destructive, non-invasive and non- contacting, fast and 
automatic measurement process of the water contents and 
temperature distribution of dielectric bodies requiring 
minimal human intervention. The measurement process is 
10 insensitive to changes in product size, form and positioning. 

Other features of the current invention will become more 
apparent in the following detailed description of the 
preferred embodiment which by means of example illustrate the 
principles of this invention. 

15 Brief Description of the Drawings 

Figure 1 is a schematic diagram of a first embodiment of a 
device according to this invention. 

Figure 2 is a chart indicating a dielectric model for chicken 
anticipated for reduction of the amount of unknown variables 
20 of the sample's dielectric behaviour. 

Figure 3 is a chart indicating a dielectric model for bread 
anticipated for reduction of the amount of unknown variables 
of the sample's dielectric behaviour. 

Figure 4 illustrate a cross section of a bread loaf, where the 
25 dielectric model from figure 3 is mapped. 

Figure 5 is a schematic chart indicating the evaluation 
process in order to obtain moisture, density and temperature 
data from dielectric properties. 
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Figure 6 is a schematic diagram of a second embodiment of a 
device according to this invention. 

Figure 7 is a flow chart of the whole calculation process. 

Detailed Description of Preferred Embodiments 

As shown in Figure 1 the primary elements of a measurement 
device 10 according to a first embodiment of the present 
invention are a microwave generator 11, a transmitting antenna 

12, a receiving antenna 13, an analyser 14. These elements 
work together to analyse the distribution of material 
properties (such as water contents, density and temperature) 
in a material sample 16. The sample is carried on a conveyor 
means 17, which may consist of a slide table mounted on a 
linear motor, and is arranged in a measurement gap between 
said transmitting antenna 12 and receiving antenna 13 . 

The generator 11 is connected to the transmitting antenna 12 
and generates electromagnetic radiation, which is transmitted 
from the transmitting antenna 12 towards the receiving antenna 

13, The material sample 16 is placed between said transmitting 
antenna 12 and said receiving antenna 13, which indicate that 
at least a part of the transmitted radiation passed through 
the material sample 16. The electromagnetic radiation is 
transmitted in the form of signals 18, each having a first 
amplitude and phase, and a different frequency within a 
frequency range. 

The generator 11 is also connected to the analyser 14, and 
information regarding the amplitude and frequency of each 
transmitted signal 18 is sent to the analyser 14. 

The transmitted signals 18 pass, at least partially, through 
the material sample 16 and are received by the receiving 
antenna 13 as receiving signals 19 each having a second 
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amplitude and phase, which may be different from the first 
amplitude and phase, for each different frequency. 

The receiving antenna 13 is connected to the analyser 14, 
which receives information regarding the received signals 19. 
The analyser 14 compares the amplitude and phase of the 
transmitted signal with the corresponding amplitude and phase 
for the received signal, for each transmitted frequency. 

Each transmitting antenna 12 is designed to emit 
electromagnetic radiation of a set of selected frequencies 
partially impinging on and flowing through the material 
samples 16. Each receiving antenna 13 is designed to receive 
electromagnetic radiation emitted from any transmit antenna 12 
and at least partially transmitted and reflected by the 
material sample 16. The receiving antenna 13 may be set up at 
one or more positions enabling to scan the material sample 16. 

The analyser 14 acts as interface between the raw data and the 
user. The output of the analyser 14 consists of a three- 
dimensional picture of the material sample's properties as 
density, water contents and/or temperature. 

Information about the microwave attenuation and runtime (or 
phase and damping of the microwave power wave) between the 
transmitting antennas 12 and receiving antennas 13 are 
calculated in the analyser 14. For each frequency of the 
chosen frequency set and for a chosen set of transmitting- 
/receiving antenna pair and at a fixed point on the material 
sample 16 such a calculation is performed. 

In this embodiment of the invention it is assumed that the 
shape of the material sample 16 is known, and a three 
dimensional image of the material sample is stored in a memory 
15 connected to the analyser 14 . The three dimensional image 



may be used to calculate cross -sectional images for each 
measurement position of the material sample on the conveyor 
means 17. Examples of a material where the three dimensional 
image is known are fluids passing through the gap in a tube or 
samples having a defined shape, such as candy bars. 

For all measurement positions along the material sample 16, 
the results of the damping and phase measurement, for all 
frequencies, are used to determine an electromagnetic picture, 
which is obvious for a person skilled in the art, and since 
this is not an essential part of the invention these steps are 
not disclosed in this application. The position information 
from the memory is saved as a three dimensional surface 
position data set describing the three dimensional contour of 
the material sample 16. 

The material properties (such as water contents, density and 
temperature) in a material may be obtained by interpolation of 
the material property distributions in the following. 

Assume a set of material samples has been measured previously 
as references. The data sets are stored in their original size 
or in a transformed form to reduce the data size. For these 
materials, the distribution of the parameters to be measured 
is known. These can be different temperatures, different 
temperature profiles, different density and water contents 
distributions. Extracted parameters of the measurement of 
these reference products form a point in a high dimensional 
vector space. To each point in this space a specific 
distribution of the parameters to be determined is associated 
by interpolation of the adjacent points of the reference 
measurements. The measurement results on an unknown product is 
now associated with another point on this vector space. Since 
the parameter distribution to be measured is known for a 
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certain region in the vector space, the distribution 
associated with the measured point yields the measurement 
result . 

On the other hand direct calculation of the material property 
distribution may be applied. 

Together with a three dimensional model of the dielectric 
structure of the material sample this three dimensional 
picture is used to determine regions within the measurement 
gap where the (yet unknown) dielectric function of the 
material can be assumed non-changing. Figure 2 illustrates a 
model for chicken 2 0 and figure 3 illustrates a model for 
bread 3 0 . 

Each model comprises several regions 21, 31, where the 
dielectric function is assumed to be constant. The number of 
regions in the models may be adjusted, even during the process 
of obtaining the material properties, to obtain a smooth, but 
not too smooth, curve for the dielectric constant as a 
function of x and y co-ordinates, e(x,y) . 

The regions in figure 3 are divided by concentric circles 32 
and a number of mapping points P1-P14 are arranged on the 
outer concentric circle 33. The distance between each mapping 
point is preferably essentially equal . 

The appropriate model is adapted to the three dimensional 
image of the sample material, in this example a bread loaf. 
Figure 4 illustrates a cross-section 40 of a three dimensional 
image of the bread loaf together with an x-axis and an y-axis. 
The contour of the bread is indicated by the line 41, which is 
derived from the three dimensional surface position data set 
stored in the memory, and the mapping points P1-P14 in figure 
3 are mapped upon the contour line 41. The concentric circles 
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32 in figure 3 are adjusted after the shape of the contour 
which is illustrated by the lines 42 in figure 4 divides the 
cross section of the bread loaf into regions 4 3 where the 
dielectric constant is assumed constant. 

Below is described a simplified approach of CSI, anticipating 
regions where the dielectric function is constant, as 
indicated in the models described in figure 2 and 3 . 

Starting with the relation between the scattered field at a 
given location as a function of the contrast source one can 
simplify the solution process considerably when the location 
of regions where the dielectric function is constant are known 
a priori : 

uj(p)= uf (p)+k 2 \G{p iq ) x(g) uj(g) Ms) 

D 

where G denotes again the two-dimensional Green's function of 
the electromagnetic problem 

G(p, g )^HM(k .\p-g\) 

and the polarisability^ depends on the dielectric function of 
the material s being constant on the region £>,and the 
background £ b in the following way: 
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Obviously the above step reduce the matrix size from the 
number of contrast sources to the number of different regions 
taken into account. 

From the above a similar integral equation for the scattered 
5 electric field at any point r is set up. 

Fj<r)-{**±Z. • J«. W fr-k-«()-hW + «y fc, W]-*fe) 

For this relation a similar solution process as in the general 
case is applied: 

10 Below is described a calculation of the dielectric function 

for one pair of antennas for various frequencies for frequency 
independent polarisation. 

Starting with the relation between the scattered field at a 
given location as a function of the contrast source one can 
15 simplify the solution process considerably when the location 

of regions where the dielectric function is constant are known 
a priori : 

/)=«** W) + * 2 \G(p,q,f\ X {q)- u(q,f) dAa) 

D 

2 0 In a step similar to the above procedure, the relation is 
simplified by introducing regions where the dielectric 
function is assumed to be constant: 

«7(p./)-« 4>e (p./)+* a Z^ • \G{p^f) u{ q> f) dy{ q ) 

25 where G denotes again the two-dimensional Green's function of 
the electromagnetic problem 
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and the polarisability % u depends on the dielectric function 
of the material € being constant on the region D n and the 
background s b in the following way: 



From the above a similar frequency dependent integral equation 
for the scattered electric field at any point r is set up. 



For this relation a similar solution process as in the general 
case is applied. 

Below is described a calculation of the dielectric function 
for one pair of antennas for various frequencies for frequency 
dependent polarisation. 

A first order approximation for the frequency dependence of 
the polarisation is obtained by grouping the measurement 
frequencies in two groups, a group at lower and a group at 
higher frequencies. The above summarised calculation process 
is repeated twice and the difference in the obtained 
polarisation values gives a measure for its frequency 
dependence . 




The wave vector k is defined to be the wave propagation 
constant in the background medium given by s rb9 fi rb : 




In order to calculate the material parameters based on 
dielectric data, the relation between the material parameters 
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as density, temperature and water content is needed. For most 
applications the following model for the temperature 
dependence of the dielectric function of water ( extracted 
from experimental data published in IEEE Press 1995 by A. 
Kraszewski, with the title "Microwave Aquametry") is: 

An approach (based on a simple volumetric mixing relation 
yields the dielectric chart depicted in Figure 5 where the 
real and imaginary parts of the dielectric function are taken 
as independent co-ordinates : 

*( T > C «20><t) = ( l ~ CX20)' Ghosts 'd + CH20 ' fen* fr) ~ ' d) ( 2 > 

Obviously every point in the complex dielectric plane stands 
for a unique water contents and material temperature when the 
dielectric properties of the dried base material do not change 
considerably. An unique density -temperature plot is obtained, 
when the water contents is uniform. 

From the spatial distribution of the dielectric function of 
the material sample 16, its density distribution moisture 
content and temperature are readily obtained applying a water 
model (see equation 1) and a mixing relation (see equation 2) . 
This part of the evaluation is shown schematically in Figure 
5, a schematic view of the complete calculation process is 
given in Figure 7 . 

The imaginary part of the dielectric constant Im(e) forms a 
first axis in figure 5 and the real part of the dielectric 
constant Re (e) forms a second axis, perpendicular to the first 
axis. The real part is positive and the imaginary part is 
negative. Any material without water content have a specific 
dielectric constant, so called e dr y/ which vary between point 50 
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and 51 depending on the material, both only having a real 
part. On the other hand, pure water having a temperature of 4°C 
has a dielectric constant 52 comprising both a real part and 
an imaginary part, and when the temperature of the water 
5 increase it follows a curve 53 to a point where pure water has 
a temperature of 99°C and a dielectric constant 54. The real 
part of the dielectric constant for materials containing any 
amount of water decreases with higher temperature and the 
imaginary part of the dielectric constant for materials 
10 containing any amount of water increases with higher 

temperature- For illustration see the dashed lines in figure 5 
for water content of 25, 50 and 75%. 

An example of a dielectric value 55 is indicated in figure 5. 
The value 55 is situated within a region 56 delimited by the 
15 curve 53, stretching between point 52 and 54, a straight line 
between point 54 and e dry and a straight line between edry and 
point 52. As mentioned before, if the temperature increase, 
with constant water content, the value of the dielectric 
constant 55 moves to the left in the graph as indicated by the 

2 0 arrow 56, and if the temperature decrease, with constant water 

content, the value 55 moves to the right as indicated by the 
arrow 57. On the other hand, if the water content decrease, 
with constant temperature, the value 55 moves towards e dry as 
indicated by the arrow 58, and if the water content increase, 
25 with constant temperature, the value 55 moves away from e dry as 
indicated by the arrow 59. 

For each defined region 43 the calculated, or estimated, 
dielectric constant may be directly transformed into water 
content and temperature. 

3 0 Figure 6 illustrates a measurement device 60 according to a 

second embodiment of the present invention. This embodiment 
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comprises the same parts as the first embodiment described in 
connection with figure 1, except that the memory 15 is 
replaced with a video imaging arrangement comprising two video 
cameras 61 and 62, both connected to an evaluation unit 63 , 
5 which in turn is connected to the analyser 14. 

Each video camera 61, 62 continuously take pictures of the 
material sample 16. The pictures are sent to the evaluation 
unit 63, where a three dimensional picture is created using 
known techniques. The resulting three dimensional picture 
10 similar to the one that was stored in the memory 15 in the 
first embodiment . 

By using video imaging the system gets more flexible and it is 
possible to use the measurement device on material samples 
having an unknown shape or even a changing shape depending the 
15 water content and/or the temperature. 

In the above evaluation the major reason to use video imaging 
is to reduce the number of unknowns in the calculation process 
to obtain the dielectric function's distribution in the 
material sample. The obtained reduction in calculation time is 

20 necessary (at least in today's available calculation power) to 
speed up the measurement process . In this preferred 
embodiment, the material samples are easily accessible to 
video imaging. If this is not the case, alternative solutions 
are ultrasound imaging. If the material samples have a simple 

25 geometric form or if subsequent material samples are very 

similar, no extra imaging is necessary to perform the above 
calculation process as described in the figure 1. 

The calculation of the dielectric image (of a two-dimensional 
cross section ) of the material sample in the measurement gap 
3 0 is accomplished by solving the previously described inverse 
scattering problem. 
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Both video cameras 61 and 62 image the part of the measurement 
gap. The location of the cameras 61 and 62 are chosen in a way 
to enable the reconstruction of a three-dimensional picture 
where the material sample 16 is positioned within the 
measurement gap . 

For each position of the material sample 16 a three- 
dimensional picture of the sample location in the measurement 
gap is calculated based on images taken by the video cameras 
61 and 62. 

In addition the position information contained in the optical 
image is used together with a priori knowledge of the material 
structure the obtain a first guess of the dielectric structure 
under measurement. This enables to reduce the number of 
unknowns of the dielectric imaging calculation process 
drastically (about two orders of magnitude) and to speed up 
the calculation considerably. 

Figure 7 show a schematic view of the complete calculation 
process for the device according to the invention. 

As previously described in connection with figure 1 and figure 
6, the input data to the analyser comprises the microwave 
transmission measurements, i.e. information regarding the 
emitted signals 18 (amplitude and phase for each used 
frequency) and the detected signals 19 (amplitude and phase 
for the corresponding frequency) . This information is input in 
the calculation process, 71. 

Information regarding the image contour of the material sample 
16 is also needed and inputted into the process, 72. A 
predetermined resolution of the image contour is used to start 
the calculation process. The resolution may be increased or 
decreased dependent on the calculation results, as described 
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below. Information regarding the position of the material 
sample 16 in the measurement gap is also inputted into the 
process in 72 . 

The information from 72 is used to establish an object 
5 geometry, 73. A model, for instance as described in figure 3, 
is thereafter used to determine regions wherein the dielectric 
function is assumed of the first order, i.e. constant. The 
number of regions used is set in the model. The selected 
model, in this case model 30, is used to establish regions in 
10 the material sample 16 by adjusting the concentric circles to 
the result of the object geometry from 73, which is done in 
74, as described in figure 4. 

The geometry assumptions from 74 and the result from the 
microwave transmission measurement from 71, is thereafter used 
15 to calculate the dielectric constant within each region, 75. 

The calculation process have previously been described in this 
application. 

Another piece of information is needed to convert the 
dielectric constant into water content and/or temperature, 
2 0 that is the dielectric constant for the material sample 16, 
when there is no water content in the material, E^ry. This 
information may be obtained from literature or from previously 
made measurements on similar material samples, 76. 

This information is used to establish the equations defining 
25 the relation between the dielectric constant and the water 
content and temperature, as described in connection with 
figure 5, 77. 
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The resulting dielectric constant within each region from 75 
is thereafter translated (or converted) into water content and 
temperature , 78. 
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If the calculation process is well established the following 
steps may be unnecessary, but in most cases they are necessary 
to avoid unreasonable results. 

In 79, a check is made to determine if the resulting 
5 temperature and water contents are reasonable, i.e. the 

temperature is greater than zero, T>0, the water content is 
greater than zero, C H2 o>0 (i.e. Im(e)<0) and if the water 
content is less than 100 %, Ch2o<100%. 

If any of the above mentioned checks does not pass, the 
10 calculation process is fed back via 80, where the position of 
the material sample is updated. If video cameras are used, as 
described in figure 6, a new image contour of the material 
sample is used to repeat the steps 74, 75 and 78. In the case 
where the image contour information is previously stored in a 
15 memory, as described in figure 1, the calculation process may 
make a small adjustment to the material size, deform the 
material contour, translate the material in one direction and 
the repeat steps 74, 75 and 78. 

If no objections are raised regarding reasonable results in 79 
20 the process continue to 81, where the smoothness of the curve 
describing the dielectric function across the cross section of 
the material sample is investigated. If the dielectric 
function is too smooth or not enough smooth, the process is 
fed back via 82, where the resolution of the image contour is 
25 changed. Thereafter the steps 73, 74, 75 and 78 are repeated 
before the checks 79 and 81 are performed again. 



There is also a possibility to change the number of regions in 
the used model in 74 to increase or decrease the number of 
regions to calculate. 
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If the smoothness of the curve is acceptable, the process 
proceed to 83, where a new dry dielectric constant, e dry , of the 
material is calculated depending on the calculated results in 
the process. If the calculated dry dielectric constant, edry/ 
5 does not correspond with the used dry dielectric constant, 
£dry. priori the dielectric constant is updated in 84 and the 
translation of the dielectric function in step 76, 77 and 79 
are repeated, before the checks 79, 81 and 83 are performed 
again* 

10 If no objections are raised in 83, the process present the 
results in the form of water content and or temperature. 

The calculation process described in figure 7 is normally 
performed for a position of the material sample in the 
measurement gap. When the calculation process is completed the 
15 conveyor means 17, on which the material sample 16 is moved to 
a new position where another measurement is performed. The 
updated inf ormat ion, regarding , Edry/ number of regions, 
position of material, and so on, are used at the next position 
to speed up the process . 

2 0 In a further embodiment of the present invention a multiple of 
receiving antennas may be used to allow a single processing , 
as described in figure 7, to establish the three dimensional 
temperature , or water content, distribution within the 
material sample. 

25 The calculation process in figure 7 only describe the 

embodiment where the model is used to establish regions, where 
the dielectric constant is assumed constant- 
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Claims 

1. A device (10; 60) for measuring the distribution of 
selected properties of a material (16) arranged on a conveyor 
means (17) , said device comprises 

5 an emitter (12) of electromagnetic radiation arranged 

at one side of said conveyor means (17) , said emitter emits 
electromagnetic radiation (18) in a multiple of frequencies in 
a selected frequency range towards said material (16) , 

at least a sensor (13) arranged on an opposite side 
10 of said conveyor means (16) , compared to said emitter (12) , 
said sensor (13) detects electromagnetic radiation (19) in 
said selected frequency range being emitted from said material 
(16) , said electromagnetic radiation originating from said 
emitter (12) , and 

15 - an analyser (14) is arranged to receive information 

regarding said emitted electromagnetic radiation (18) and said 
detected electromagnetic radiation (19) , said information 
comprises amplitude and/ or phase for each selected frequency, 

characterised in that said device further comprises: 

20 - an image device (15; 61-63) which is arranged to send 

information to said analyser (14), to create a three 
dimensional contour of the material, 

said analyser (14) is arranged calculate the selected property 
distribution in said material (16) based on the received 
2 5 information . 

2. The device according to claim 1, wherein said image 
device is a memory (15) , having stored information regarding 
the three dimensional contour of the material (16) . 
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3. The device according to claim 1, wherein said image 
device comprises at least an imaging sensor (61, 62) is 
connected to an image processing device (63) , said at least 
one imaging sensor (61, 62) each detects an image of said 

5 material, which is processed into a three dimensional contour 
of said material in said image processing device (63) . 

4 . The device according to claim 3 , wherein said at 
least one imaging sensor (61, 62) detects a picture of the 
reflectivity in optical wavelengths. 

10 5. The device according to claim 4, wherein said at 

least one sensor of a second type is a video camera (61, 62) . 

6. The device according to any of claim 3, wherein said 

imaging sensor detects a picture of the reflectivity and 
transmissivity and propagation speed of sound waves. 

15 7. The device according to claim 6, wherein said at 

least one sensor of a second type is an ultrasound imaging 
device. 

8. The device according to any of claims 1-7, wherein 
said analyser (14) is provided with means to interpolate 

2 0 previously measured results, stored in a memory, to obtain the 
selected property distribution in said material (16) . 

9. The device according to any of claims 1-7, wherein 
said analyser (14) is provided with means to calculate the 
dielectric distribution in said material (16) and convert said 

25 dielectric distribution into the selected property 
distribution in said material (16) . 



10. The device according to claim 9, wherein said means 

to calculate the dielectric distribution comprises 
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a three dimensional model (20 , 30) determining 
regions within said material (16) where the dielectric 
function is assumed non- changing, and 

means to apply said model (20, 30) to said three 
5 dimensional contour of the material, 

whereby a dielectric distribution is obtained. 

11. The device according to claims 9 or 10, wherein said 

device is provided with means to convert said dielectric 
distribution into the selected property distribution. 

10 12. A method for measuring selected properties of a material 
(16) arranged on a conveyor means (17), said method comprising 
the steps: 

(a) emitting electromagnetic radiation (18) in a selected 
frequency range towards said material (16) from an emitter 

15 (12) arranged at one side of said conveyor means (17), 

(b) detecting electromagnetic radiation (19) in said selected 
frequency range in a sensor (13) arranged on an opposite side 
of said conveyor means (16) , compared to said emitter (12) , 
said electromagnetic radiation (19) being emitted from said 

2 0 material (16) , said electromagnetic radiation (19) originating 
from said emitter (12) , 

(c) transmitting information, comprising amplitude and/or 
phase for each selected frequency, regarding said emitted 
electromagnetic radiation (18) and said detected 

25 electromagnetic radiation (19) to an analyser (14) , 



characterised in that said method further comprises the steps: 
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(d) transmitting information, comprising information to 
create a three dimensional contour of the material, to said 
analyser (14) , and 

(e) calculating the selected property distribution in said 
5 material (16) by analysing said information in the analyser. 

13. The method according to claim 12, wherein the step of 
calculating the selected property distribution comprises the 
step of interpolating previously measured results, stored in a 
memory . 

10 14. The method according to claim 12, wherein said step of 
calculating the selected property distribution comprises the 
steps : 

calculating of the dielectric distribution in said 
material using the information regarding emitted (18) and 
15 detected (19) electromagnetic radiation, and 

converting the dielectric distribution into the selected 
property distribution in said material (16) . 

15. The method according to claim 14, wherein said step of 
calculating the dielectric distribution comprises the steps: 

2 0 - providing a three dimensional model (20, 30) determining 

regions within said material (16) where the dielectric 
function is assumed non-changing, and 

applying said model (20, 30) to a three dimensional 
contour obtained by an imaging device (15; 61-63) . 

25 16. A system comprising a device according to any of claims 1 
to 11, said instrument evaluates the measurement data based on 
a method according to any of claims 12 to 15 in order to 
obtain information on the spatial distribution of the said 




dielectric and magnetic properties of said materials, which 
information is used to calculate the temperature, density 
and/or water contents distribution of said materials. 
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Summary 

The present invention relates to a device for measuring the 
distribution of selected properties of materials, said device 
comprises an emitter of electromagnetic radiation and 
5 furthermore at least one sensor of a first type. The emitter 
emits electromagnetic radiation in a selected frequency range 
towards said materials and a sensor of the first type detects 
electromagnetic radiation in a selected frequency range coming 
from said materials* The detected electromagnetic radiation 

10 having been emitted by said emitter. The device also comprises 
means to generate a three dimensional image contour 
information regarding the said material's position in space, 
and an analyser which (a) receives information from said 
sensors and (b) processes this information and (c) generates 

15 signals containing information about the distribution of said 
properties as output. The invention also relates to a system 
and a method for measuring the distribution of selected 
properties of materials. 
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